Getting the functions and the data
library(dispRity)
source("functions.R") ; load.functions(test = FALSE)
## NTS.R :
## apply.test.R :
## char.diff.R :
## cor.matrix.R :
## export.ntw.R :
## get.aov.data.R :
## modify.matrix.R :
## plot.char.diff.density.R :
## plot.cor.matrix.R :
## plot.results.R :
## read.TreeCmp.R :
## write.nexus.std.R :
dyn.load("../Functions/char.diff.so")
load("../Data/NTS/matrix.nts.Rda")
source("../Data/load.data.R")
Comparison to ‘best’ tree (RF normalised distance)
op <- par(mfrow=c(3,1), mar = c(2, 4, 1, 2))
plot.results.single(t25_list, metric = 2, what = "best",
legend = TRUE, NTS = TRUE, xlab = "", ylab = "RF for 25 taxa",
axislab = FALSE)
plot.results.single(t75_list, metric = 2, what = "best", ylab = "RF for 75 taxa",
xlab = "", axislab = FALSE, NTS = TRUE)
plot.results.single(t150_list, metric = 2, what = "best",
ylab = "RF for 150 taxa", NTS = TRUE)
par(op)
Comparison to ‘best’ tree (Triplets normalised distance)
op <- par(mfrow=c(3,1), mar = c(2, 4, 1, 2))
plot.results.single(t25_list, metric = 4, what = "best",
legend = TRUE, NTS = TRUE, xlab = "", ylab = "Tr for 25 taxa",
axislab = FALSE)
plot.results.single(t75_list, metric = 4, what = "best", ylab = "Tr for 75 taxa",
xlab = "", axislab = FALSE, NTS = TRUE)
plot.results.single(t150_list, metric = 4, what = "best",
ylab = "Tr for 150 taxa", NTS = TRUE)
par(op)
Comparison to ‘null’ tree (RF normalised distance)
op <- par(mfrow=c(3,1), mar = c(2, 4, 1, 2))
plot.results.single(t25_list, metric = 2, what = "null",
legend = TRUE, NTS = TRUE, xlab = "", ylab = "RF for 25 taxa",
axislab = FALSE)
plot.results.single(t75_list, metric = 2, what = "null", ylab = "RF for 75 taxa",
xlab = "", axislab = FALSE, NTS = TRUE)
plot.results.single(t150_list, metric = 2, what = "null",
ylab = "RF for 150 taxa", NTS = TRUE)
par(op)
Comparison to ‘null’ tree (Triplets normalised distance)
op <- par(mfrow=c(3,1), mar = c(2, 4, 1, 2))
plot.results.single(t25_list, metric = 4, what = "null",
legend = TRUE, NTS = TRUE, xlab = "", ylab = "Tr for 25 taxa",
axislab = FALSE)
plot.results.single(t75_list, metric = 4, what = "null", ylab = "Tr for 75 taxa",
xlab = "", axislab = FALSE, NTS = TRUE)
plot.results.single(t150_list, metric = 4, what = "null",
ylab = "Tr for 150 taxa", NTS = TRUE)
par(op)
Plotting the results for 25 taxa
## Comparisons to the best tree (RF)
plot.results.single(t25_list, metric = 2, what = "best",
main = "comparison to 'best' tree (25 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the best tree (Triplets)
plot.results.single(t25_list, metric = 4, what = "best",
main = "comparison to 'best' tree (25 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the null tree (RF)
plot.results.single(t25_list, metric = 2, what = "null",
main = "comparison to 'null' tree (25 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the null tree (Triplets)
plot.results.single(t25_list, metric = 4, what = "null",
main = "comparison to 'null' tree (25 taxa)", legend = TRUE, NTS = TRUE)
Plotting the results for 75 taxa
## Comparisons to the best tree (RF)
plot.results.single(t75_list, metric = 2, what = "best",
main = "comparison to 'best' tree (75 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the best tree (Triplets)
plot.results.single(t75_list, metric = 4, what = "best",
main = "comparison to 'best' tree (75 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the null tree (RF)
plot.results.single(t75_list, metric = 2, what = "null",
main = "comparison to 'null' tree (75 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the null tree (Triplets)
plot.results.single(t75_list, metric = 4, what = "null",
main = "comparison to 'null' tree (75 taxa)", legend = TRUE, NTS = TRUE)
Plotting the results for 150 taxa
## Comparisons to the best tree (RF)
plot.results.single(t150_list, metric = 2, what = "best",
main = "comparison to 'best' tree (150 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the best tree (Triplets)
plot.results.single(t150_list, metric = 4, what = "best",
main = "comparison to 'best' tree (150 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the null tree (RF)
plot.results.single(t150_list, metric = 2, what = "null",
main = "comparison to 'null' tree (150 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the null tree (Triplets)
plot.results.single(t150_list, metric = 4, what = "null",
main = "comparison to 'null' tree (150 taxa)", legend = TRUE, NTS = TRUE)
Checking the results for 25 taxa
## Comparisons to the true tree (RF)
plot.results.single(t25_list, metric = 2, what = "true",
main = "comparison to 'true' tree (25 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the true tree (Triplets)
plot.results.single(t25_list, metric = 4, what = "true",
main = "comparison to 'true' tree (25 taxa)", legend = TRUE, NTS = TRUE)
## Checking comparisons (RF)
plot.results.single(t25_list, metric = 2, what = "check",
main = "Check comparison (25 taxa)", legend = TRUE, NTS = TRUE)
## Checking comparisons (Triplets)
plot.results.single(t25_list, metric = 4, what = "check",
main = "Check comparison (25 taxa)", legend = TRUE, NTS = TRUE)
Checking the results for 75 taxa
## Comparisons to the true tree (RF)
plot.results.single(t75_list, metric = 2, what = "true",
main = "comparison to 'true' tree (75 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the true tree (Triplets)
plot.results.single(t75_list, metric = 4, what = "true",
main = "comparison to 'true' tree (75 taxa)", legend = TRUE, NTS = TRUE)
## Checking comparisons (RF)
plot.results.single(t75_list, metric = 2, what = "check",
main = "Check comparison (75 taxa)", legend = TRUE, NTS = TRUE)
## Checking comparisons (Triplets)
plot.results.single(t75_list, metric = 4, what = "check",
main = "Check comparison (75 taxa)", legend = TRUE, NTS = TRUE)
Checking the results for 150 taxa
## Comparisons to the true tree (RF)
plot.results.single(t150_list, metric = 2, what = "true",
main = "comparison to 'true' tree (150 taxa)", legend = TRUE, NTS = TRUE)
## Comparisons to the true tree (Triplets)
plot.results.single(t150_list, metric = 4, what = "true",
main = "comparison to 'true' tree (150 taxa)", legend = TRUE, NTS = TRUE)
## Checking comparisons (RF)
plot.results.single(t150_list, metric = 2, what = "check",
main = "Check comparison (150 taxa)", legend = TRUE, NTS = TRUE)
## Checking comparisons (Triplets)
plot.results.single(t150_list, metric = 4, what = "check",
main = "Check comparison (150 taxa)", legend = TRUE, NTS = TRUE)